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Abstract 

In this paper, a new method for generation of infinite series of symmetric 
identities written for exponential sums in real numbers is proposed. Such sys- 
tems have numerous applications in theory of numbers, chaos theory, algo- 
rithmic complexity, dynamic systems, etc. Properties of generated identities 
are studied. Relations of the introduced method for generation of symmetric 
exponential sums to the Morse-Hedlund sequence and to the theory of magic 
squares are established. 

Key Words: Real exponential sums, generation, Morse-Hedlund sequences, magic 
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1 Introduction 



Let us consider two sets of numbers, {xi : I < / < A:} and {j, : I < / < ^} such that 
the following system of identities is satisfied. 



~\~ X2 ~l~-^3 ~1~ ••• ~^ -^k—l ~^ "^k 

x^ ~\~ x^ -\- x^ ~^ -^j^—x ~^ -^ir 



x'{+x1+xl + ...+xl_^+xl 



y\+y\+y\ + -+y\-i+y\ 
y\+yl+yl + -+yl-i+yl 
y\+yl+yl + -+yl-i+yl (i) 

y'{+yl+yl + -+yU+yl 



Systems of the type ([T]) have a significant importance in mathematics. Such 
systems were widely studied in the number theory (see H [13] [l8]|23l). For exam- 
ple, the Tarry -Escott problem (see, for example, (T TtI) deals with finding disjoint 
sets {xj : \ <i<k} and {yi : I < / < ^} of integers satisfying ([T|l. The systems ([B 
have also a straightforward connection to the Hilbert-Kamke problem for integers 
(see ll8j[T0l[Tl]) being an extension of the famous Waring problem (see ||24|| ). 

Identities ([T]) considered as a description of the process of exponential sums ex- 
pansion have close relations to the chaos theory (see [ 16 1) and its applications, for 
instance, in economics (see |i5j|). Integer solutions to Q in the context of integrals 
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of piece-wise constant functions, emergent calculations, and integer code series 
can be viewed (see 1 14l |T5l) in terms of structural complexity of systems and al- 
gorithms. Particularly, studying structural complexity of multiextremal functions 
gives new tools for solving global optimization problems being a new actively rais- 
ing research area (see ||9] [14] [T5ll22l '). 

Together with the system ([T]) the following system dJ]) can be considered. 



1-1 +-^ei 



X, +XT+Xn + ... + 



3 3 



3 

Xn ~1~ . . . ~1~ Xf\ 



+ x^ 
+ x\ 



y x\+xi+xi + ...+4_^+xi 



y\+y\+y\ + ■■■+yl,-i+yl, 
y\+yl+yl + ■■■+yl^-i+yl^ 
y\+yl+yl + ■■■+yl,-i+yl, 

y'l+y2+y'l + -+yl-i+yl 



(2) 



Such systems have interesting connections (see fTSl) to algorithmic complexity 
(see ||6j[i2,.21J) and period doubling in physical systems (see Ii7,ii5j)- In this case, 
the numbers 9i , 62 , . . . , 6„- 1 , 6„ = from ([T} become 



62 — 261, 63 — 262, 6„— 26„_i 

and, for example, for 61 =2 the following piramide is obtained. 



(3) 



Jfj -\- X2 H~ -^3 "h • • • ~\~ Xj -\- X^ 



X] ~\- X9 H" -|- . . . -|- j^l 



k-\ 



y\+y2 

y\+yl+yl+yl 

y\+yl+yl + -+y^i+yl (4) 



y\+yi+y"^ + -+yU+yln 



In this paper, for any integer power n, a method for generation of two sets {xi : 
\<i<k} and {j,- : 1 < / < ^} of real numbers satisfying both ([T} and (|2]l for certain 
numbers 61 < 62 < . . . , 6„_i < 6„ = is proposed. Particularly, piramidale systems 
of the type ^ are considered. Properties of the introduced sets {xj : \ <i<k} 
and {yi : 1 <i <k} are studied. It is shown that the sets found for a given n are a 

basis for the sets satisfying ^ and (|2]) with the powers n + \,n + 2, In order 

to obtain solutions to a system having a higher power ?i + 1 it is necessary to add 
certain numbers to the sets {xf : \ <i<k} and {y,- : \ <i<k} found for the power 
n. In such a way, by increasing n, it is possible to speak about sequences {jc,} and 
{yi] satisfying a sequence of systems ([T|l and (O. 

The next section contains the main results. A few examples are given in the 
last section. It also presents relations of the introduced method for generation 
symmetric exponential sums to the Morse-Hedlund sequences (see |19|) being a 
particular case of the Prouhet sequences (see [20 J ) and to the theory of magic cubes 
developed in HI 13. 
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2 Main results 



In this section, two sets of real numbers {xt '■ I <i < fcjand {yt : I <i <k} satis- 
fying ([T]) and (|2l) for any given power n are constructeqj. Let us start by choosing 
four real numbers a,b,c, and d satisfying the following simple equality 

a + b = c + d. (5) 

Given numbers a,b,c,d from ^ and any real numbers ki,k2, ■ ■ . ,kn,n > 2, 
let us introduce by recursion functions U {k\,k2,. ■ ■ , k„^ \ ) and Rj{ki,k2,..., k„^ i ) , 
where j = n-\,n. The functions L\ki),R^ki),L^{ki),R^{ki) and {ki M),^^ (ki , ^2) 
are defined as follows: 

L\ki) = {a + ki)^ + {b + ki)^ =x\+x\, (6) 
where x\ = a + ki and x\ =b + ki, 

R\ki) = {c + kiY + {d + kiY =y\ +yi (7) 
where y\ =c + ki and y\ = d + ki, 

L^{ki) = {a + kif + {b + kif=x]+xl 

R^{k,) = {c + hf + {d + hf=y\+yl. 
Then, setting X'i = c + k\+k2,Xix = c + k\+k2, and R?{k\ +^2) =x\+x\ we define 
L^{ki,k2)=l^^{ki)+R^{h+k2)=x\+xl+xl+xl = 

{a + kif + {b + kyf + {c + ki+k2f + {d + kY+k2f. 

Analogously we introduce 3^3 = a + k\ + k2,yA = b + ki + k2, and L^{k\ + k2) = 
y\ + y\ to define 

R^{k,M)=R'{ki)+L^{ki+k2)=y\+yl+yl+yl = 

(c + kif + {d + h f + {a + h+ k2f + {b + ki+ ^2)^. 

Now, using already determined functions, we can introduce by recursion func- 
tions U{k\,k2, ■ . ■ ,kn-i) where j = n — l,n: 

2" 

L"{ki,k2,...,kn-i,kn) = L"{ki,k2,...,kn-i)+R"iki+kn,k2,...,kn-i) = Y^^i, 

(=1 

(8) 

where 

2"— 1 

L'\k,,k2:---,kn-y)=Y.x'l, (9) 
(=1 



'An alternative description of thiese sets proposed by the unknown referee and is given in the 
Appendix. 
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and the numbers Xi,l < i < 2"^ , have akeady been determined previously for 
powers 1, ... ,n — 1 and the remaining numbers Xi,2"^^ < i < 2", are calculated as 
follows 

Xi = yj + k,,, i = 2"-^ +j, l<i< 2"-\ (10) 
The functions Rj{ki,k2, . ■ .,kn-\),j = n — l,n, are defined by a complete analogy: 

2" 

R'\ki,k2,...,k„^i,li„) = R"{ki,k2, . . . ,k„^i) + L"{ki + kn,k2, . . . ,k„^i) = Y.y'^ , 

i=l 

(11) 

R"{kuk2,...,k„-i)=Y.yi, (12) 

i=l 

where j, , 1 < / < 2"^\ have already been determined and 

yi = Xj + kn, i = r-'+j, l<i<2"-'. (13) 

Let us write for illustration the functions L^{ki,k2,kj) and R^{ki,k2,kT,) in their 
explicit form. 

L\ky,k2M)=L\k^,k2)+R\h +k2M) = 

= {a + hf + {b + hf + {c + ki+k2f + {d + k\ +^2)^ + 
+{c + ky+kif + {d + ki+krif + {a + ki+k2 + hf + {b + ky+k2 + hf, 
R^{ki,k2,k3)=R\ki,k2)+L^{ki+k2,k3) = 
= (c + ki)^ + {d + ki)^ + {a + ki+k2)^ + {b + ki+k2f + 
+{a + ki+k3f + {b + ki+k3f + {c + ki+k2 + k3)^ + {d + ki+k2 + k3)\ 
Now we are ready to formulate the first theorem. 

Theorem 1 If L\0) = R^{0), i.e., (O holds, then for all n>\ and for all real 
numbers ki,k2, ■ ■ ■ ,k„ it follows 

L"{kuk2,...,kn)=R'\k,,k2,...,k„). (14) 

Proof. Theorem is proved by induction. For n = 1 we obtain from Q, and O 
thatLi(A:i) =R^iki). 

Suppose now that (fT4l) holds for / such that 1 < / < « — 1 and prove that this 
assumption implies truth of ([T4l) for / = n. To proceed we need some designations. 
We denote by card{^^ k^} the number of items A:^, 1 < m < n, in the sum Y^m^m- 
We also introduce sets 5,y, 1 < / < ?i,0 < 7 < /, as follows 

Sio = {0}, 

"^'7 = {■^'7 • ^ij = Y^k,n,card{Y^k,n} = j,n-i <m<n}, 0<j< i. (15) 
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Thus, a set Sij contains as elements all possible sums Y,m^m,n — i <m <n, having 
j items. 

Now we are ready to calculate the left part of (fT4] |. i.e., L''{ki,k2, . ■ ■ ,kn). By 
using ^ and (flOl) we obtain 

L" {k.M,..., kn^i, k,)=L"{ki,k2,---,kn^i) + Y. (yj + ^«)"- (16) 

7=1 

By using Newton's binomial formula for every item {yj + k„)",l < j < 2"^^, in 
([T6l ) and then applying ([T2l) and designations ([15] ) we have 

L"{ki,k2,...,k„^i,kn) = L"{ki,k2,. ■ . ,kn-i)+R"iki,k2,. . ■ ,kn-i)+ 

"t'^ni"' ) it I /?"-'('tl+^,7,^2,^3,...,-tn-0]+2"-X", (17) 

(=1 V / j=o SijeSij 



where 




are binomial coefficients. 



We calculate R" {k\,k2, ■ ■ ■ ,kn-\-,kn) by a complete analogy. It follows from 
(E]) and ^ that 

R\k,M,---,kn-ukn)=R''{kuk2,...,kn-i)+Y,{xj + k„)\ (18) 

;=i 

Applying again in ([T8] ) Newton's binomial formula for items (xy + ^„)", 1 < j < 
2"^' , and then using ^ and (fT5l) we can write 

R''{kuk2,---,K-uK)=R'\ky,k2,---,kn-\)+L"{ky,k2,---,K-\) + 

(")(t I L"-'{ky+Sij,k2M,.--,kn-i)]+r-'K, (19) 
i=l V W ;=0 SijeSij 

Since numbers n — i we. such that 1 < n — i < n then, due to our assumption, the 
following identities 

L"^' {ki + ^2 , ^3 , • • • , = (^1 + Sij,k2 , ^3 , • • • , 

hold for all < j < / — 1 , 1 </<« — !. This result considered together with ([TtT i 
and (fT9l ) completes induction and proves the theorem. ■ 

Now we are ready to obtain the main result concerning systems of the type O 
for the case k = 2". 

Theorem! If L^{0) =R^{0), i.e., (15]) holds, then for all n>\ and for all real 
numbers ki,k2, ■ ■ ■ ,kn the numbers Xj,yj, 1 < j <2" , defined by (16]), ([71), f liOl ). and 
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in~3l satisfy d?]) and the following system 



~\~ X2 ~\~ X'^ -|~ ... ~\~ X2n _Y ~\~ Xjn 
X-y ~\~ X2 ~\~ X"^ ~h . . . ~h -^2" — 1 ~^ -^2" 



^ X-y ~\~ X2 ~l~ -^3 H~ • • • H~-^2" — 1 ~l~ ■^2" 



yl+yi+yl + ■■■ +yln_i +yln 
yl+yl+yl + --+yl<-i+yln (20) 

yi+y2+3'3 + -+3'2"-l+3'2" 



Proof. Let us introduce auxiliary functions M'{ki,k2,. ■■ ,k„^i,k„),l < i < n, as 
follows 

M'{ki,k2,. ■ .,k„^i,kn) =x\+X2+x\ + ... +X2«-i +X2n,l <i<n. (21) 
From (03] ) and ([TT] ) we can write 

n—i 

M'{ki,k2,...,k„^i,k„) = Y, 52 L'{ki+Sij,k2,k3,...,ki). 

j=o s.jes.j 

Due to (fT4l) and (fTTI) . we can rewrite this equality as 

«-! ^ 2" 

7=0 s.jES.j ;=1 

where 1 < / < «. Thus, the theorem has been completely proved. ■ 

Let us now establish some properties of the introduced numbers Xm,ym,^ < 
m < 2", from dlOll. 

Property 1 If there exists a number i > 1 such that ki = then it follows for all 
n > i that for every number 1 <m<2", there exists such a number yp that 

— yp- 

Proof. Truth of the fact follows immediately from definitions ([TOl i and (fT3] ). ■ 

The next property establishes presence of symmetry not only with respect to 
the sign '=' but also for inner subsumes of system ([T}. 

Property 2 For subsumes of the numbers Xm,ymi 

\<m< 2", from ^ the follow- 
ing identities hold 

(/>+l)2' (p+l)2' 

L 4= L '4' 0</'<2"-'> 

j=p2i+l i=p2' + l 



Proof. This result is a straightforward consequence of theorems [T]and[2l ■ 
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Property 3 If the number a + b + 2k[ + Y!i=2 ^' ^'^ integer then 

M\ki,k2, ■ . ■ ,kn-l,kn)j^g^2 = 0- 



Proof. This result is obtained from ([6) - ([T3T l. since the number {ki,k2,-- - , kn- i,kn) 
can be written in the form 

2" n 

M\kuk2.---^K-ukn) = Y.^i = ^"'\a + b + lki + Y,h)- ■ 

(=1 1=2 

Property 4 For a fast calculation of the numbers M'{ki,k2, ■ ■ ■ ,^„_i,^„), 1 <i<n, 
from f l2iD the following recursive formula can be used 

M'{ki,k2,. . . ,kn-ukn) = 2M'{ki,k2,. . . 

" )M"-'(^i,;t2,...,^„_i)+2"-i;t;:. (22) 



Proof. Property is a straightforward corollary of ([19) and dlTl . ■ 

Theorem 2 has been proved for numbers k from ([Til equal to 2". However, it is 
possible to generate solutions to ([T|) for / 2" by a suitable choice of the numbers 
a,b,c,d,ki,k2, . ■ . , k„ in such a way that some numbers xj = j,- = and so can be 
excluded from consideration in ([T|l. It is also possible to generate 'asymmetric' 
sums in sense that the number of items not equal to zero laying on the left from '=' 
does not coincide with that number for the sum laying on the right (see examples 
in the next section). 

The following theorem extends results established above from the case of two 
items in the basic identity ([5]l to the case with items. 

Theorems If given parameters ki,k2, ■ ■ ■ ,k„ in the recursive process (tS)) - ([I3\l 
the basic identity (O is substituted by the identity 

N N 

L «i = L (23) 

then, the numbers Xi,yi generated by the scheme will be solutions to systems ([7]), 
0fork = 2"-^N. 

Proof. The theorem is proved by a complete analogy with the basic case N = 2. m 
Properties similar to the properties established for the case ([5]l can be proved 
for ^ too. 
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3 Examples and relations to the Morse-Hedlund sequence 



Let us consider the first example for tlie case « = 3 witli tlie numbers 
a = —l, b = 2.l, c = 3.4, d = —2.3, 

ki=0, k2 = \, k3 = -0.5. (24) 

Application of the introduced method to the data (l24l) gives the following piramide 
dH) of identities. 

(-l)'+2.l' = 3.4'+(-2.3)', 
(-1)2+2.12+4.42 3)2 ^ 3.42+(-2.3)2+02+3.l2, 

(-l)3+2.l3+4.43+(-1.3)3+2.93+(-2.8)3+(-0.5)3+2.63 = 3.43+(-2.3)3+o3+3.l3+(-1.5)3+1.63+3.93+(-l.l 

The sums for the first, second and third identities are equal to 1.1, 26.46, and 
111.136, respectively. 

Let us now use the data from (l24l) and construct the corresponding systems of 
the type ([l). This operation leads to two groups of identities presented below. The 
first group consists of identities for « = 2 

(-l)'+2.l'+4.4'+(-1.3)' = 3.4'+(-2.3)'+0'+3.l'. 
(-l)2+2.l2+4.42+(-1.3)2 = 3.42+(-2.3)2+02+3.l2 

where sums in the first line are equal to 4.2 and in the second to 26.46. The second 
group of identities is written for « = 3. 



(-l)'+2.l'+4.4'+(-1.3)'+2.9'+(-2.8)'+(-0.5)'+2.6' = 3.4' + (-2.3)'+0'+3.l' + (-1.5)' + 1.6'+3.9' + (-1.8)i, 
(-l)2+2.l2+4.42+(-1.3)2+2.92+(-2.8)2+(-0.5)2+2.62 = 3.42+(-2.3)2+o2+3. l2+(_ 1.5)2+1.62^3 92^(_l g)2^ 
(-l)3+2.l3+4.43+(-1.3)3+2.93+(-2.8)3+(-0.5)3+2.63 = 3.4^ + {-23)^ +0^ +3.1^+{-1.5)^ + 1.6^+3.9^ + . 

The sums for the first, second, and third identities are equal to 6.4, 49.72, and 
111.136, respectively. 

The second example is in integer numbers for n = 4. In order to obtain repeated 
numbers Xi,yi and to demonstrate the nature of a 'hidden' symmetry appearance the 
following parameters 

a = l, b = 3, c = d = 2, ki=0, k2 = l, k^ = —2, ^4 = 3 (25) 

have been chosen. The pyramid from dD for these data is shown below. 

l'+3' = 2'+2', 

12+32+32+32 = 22+22+22+42, 

P+33+33+33+03 +0^+03+23 = 23+23+23+43+(-l)3+l3 + l3 + l3, 

l'*+3'*+3*+3*+0'*+0'*+0'*+2-*+5''+5-*+5^+74+2''+44+4''+4'* = 2'*+2'*+2'*+4*+(-l)*+l-*+l4+l-*+4*+6*+6*+6*+3*+3'*+3*+5'*. 

The sums for the first, second, third, and forth identities are equal, respectively, to 
4, 28, 90, and 5320. In this record all the numbers Xi,yi have been shown. Let us 
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now eliminate from the record items equal to zero and elements Xi and yj such that 
Xi = yj. In this case the symmetry is present only in an implicit form and is not 
seen from the record: 

1^+31 = 2^+21, 
12 + 32 + 32 + 32 = 22 + 22 + 22+42, 

33 + 33+33 = 23 + 23+43 + (-l)3 + l3 + l3^ 
54 + 54 + 74^44 ^ 24 + (-l)4 + l4 + l4 + 64 + 64 



Let us consider now how does system (1201 ) look like for data (I25I ). 

l'+3'+3'+3'+0i+0'+0i+2'+5'+5'+5'+7'+2'+4'+4'+4i = 2i+2i+2i+4i+(-l)i + l' + l' + l'+4i+6'+6'+6'+3'+3'+3'+5', 

1^+32+32+32+02+02+02+22+52+52+52+72+22+42+42+42 = 22+22+22+42 + (-l)2 + l2 + l2 + l2+42+62+62+62+32+32+32+52, 
l^+3'+3^+3'+03+0'+03+23 +53+53 +5^+73 +2^+43 +43+43 = 23+23+23+43 + (-l)3 + l3 + l3 + l3_|_43_|_g3_|_63_|_g3_|_33_|_33_|_33_|_53_ 
l'*+3''+3*+3''+0*+0''+0'*+2'*+5-*+5*+5-*+7*+2-*+4''+4''+4* = 2*+2''+2*+4^+(-l)*+l''+l-*+l''+4^+6*+64+6*+3''+3*+3''+5l 

The sums are equal to 48, 208, 1008, and 5320. Deleting the explicit symmetry we 
obtain the following identities where the symmetry is hidden. 

5I+5I+7I+4I = 2i + (-l)i + li + li+6' + 6\ 

52+52 + 72+42 = 22 + (-l)2 + l2 + l2 + 62 + g2^ 
53+53+73+43 ^ 23 + (-l)3 + l3 + l3+63 + 63^ 
54 + 54 + 74+44 ^ 24 + (-l)4+l4 + l4 + 64 + 64 

The corresponding sums are equal to 21, 115, 657, and 3907. 

Let us now illustrate Theorem|3]by an example having A'^ = 3 in (l23l) and shown 
in Table[T]for the power « = 4. In this example the following parameters have been 
used: 



ai = 1, a2 = 3, a3 = 7, ci = 2, C2 = 4, 03= 5, 

^1=0, ^2 = — 1) ^3 = 1-3, k4 = —2.5. 
The last line of Table [T] presents the values 

24 24 

l^yj, 1<;<4, 

i=l i=l 

showing that the numbers x;,j,, 1 < 7 < 4, are solutions to system ([T]l for n = 4 and 
A: = 24. 

Table [T] illustrates also Property |2] In fact, it is split in boxes that contain 
groups of elements of the sequences {xj} and {yj}, 1 < 7 < 4, having equal sums. 
For example, for the power 2, the box containing the numbers Xj,x\,...,x\^ ,x\2 and 
yi,y\,-,yn,yu shows that 

xj+xl + ... +x\^ +x\2 = y^i +y\ + ••• ^y\\ +3'?2- 
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Table 1 : Example for « = 4 and the starting identity 1+3 + 7 = 2 + 4 + 5 





A 


y 


X 


y 


X 


y 


X 


y 


1 
i 


1 
i 


z 


1 
i 


A 


1 
i 


Q 
O 


1 
i 


10 


z, 


J 


4. 


Q 

y 


1 ft 


77 




81 


Z>JO 


1 

J 


n 

1 


C 

J 


A Q 


2j 


'2/1 "2 




0/1 r\i 


OZJ 


A 

4 


1 
i 


U 


1 


U 


1 
i 


U 


1 
i 


U 




'X 


2 


Q 

y 


4 


27 


8 

O 


81 


16 





A 





1 A 


JO 


A/1 


ZIO 


ZjO 


1 OQA 
IzVO 


/ 




/, J 


lU,o!7 




Q'27 


lZ,10 / 


1 1 G ^QO 1 


Z / ,yo41 


8 
o 


J, J 


4. ^ 


Z,0,W" 


1 8 4.Q 


148 877 


7Q SO? 


78Q 0481 


^41 8801 


Q 

y 


0,-5 


Q T 

o,j 


jy,oy 


/;q qq 

oo,oy 


o<;r\ r\/n 


J / 1, /o / 


1 0QA1 


4 /4j,0-5zl 


1 n 


1,-5 


/,J 


i,oy 




1 QT 

Z,ly / 


Iz, 10 / 


z,oj01 


z /,Vo41 






4 3 




1 8 49 


3"^ 937 


7Q SO? 




^41 8801 

-^T^ 1 ,oo V/ 1 


12 


7,3 


5,3 


53,29 


28,09 


389,017 


148,877 


2839 8241 


789 0481 


13 


-0,5 


-1,5 


0,25 


2,25 


-0,125 


-3,375 


0,0625 


5,0625 


14 


1,5 


0,5 


2,25 


0,25 


3,375 


0,125 


5,0625 


0,0625 


15 


2,5 


4,5 


6,25 


20,25 


15,625 


91,125 


39,0625 


410,0625 


16 


-2,5 


-1,5 


6,25 


2,25 


-15,625 


-3,375 


39,0625 


5,0625 


17 


-0,5 


0,5 


0,25 


0,25 


-0,125 


0,125 


0,0625 


0,0625 


18 


3,5 


1,5 


12,25 


2,25 


42,875 


3,375 


150,0625 


5,0625 


19 


-0,2 


0,8 


0,04 


0,64 


-0,008 


0,512 


0,0016 


0,4096 


20 


1,8 


2,8 


3,24 


7,84 


5,832 


21,952 


10,4976 


61,4656 


21 


5,8 


3,8 


33,64 


14,44 


195,112 


54,872 


1131,6496 


208,5136 


22 


-0,2 


-1,2 


0,04 


1,44 


-0,008 


-1,728 


0,0016 


2,0736 


23 


1,8 


0,8 


3,24 


0,64 


5,832 


0,512 


10,4976 


0,4096 


24 


2,8 


4,8 


7,84 


23,04 


21,952 


110,592 


61,4656 


530,8416 


Sum 


61,6 


61,6 


305,08 


305,08 


1599,724 


1599,724 


9712,6972 


9712,6972 



To conclude this example let us note, that JCe = 3^2 = 4. Elimination of this couple 
of numbers from consideration gives us a solution to dl} for « = 4 and k = 22>. We 
can repeat this procedure for pairs 

■^8 =^12 = 5.3, Xi4=Ji8 = 1.5, X24= 3^20 =2.8 

arriving to a solution to (H) for « = 4 and k = 20. 



The introduced method for generation of symmetric exponential sums has in- 
teresting relations to the Morse-Hedlund sequences fT9l (being a particular case of 
the Prouhet sequences |i20J) and to the theory of magic cubes developed in |[T]|2l. 
Let us write the first items of the Morse-Hedlund sequence (see Table ^ and nu- 
merate its items by natural numbers. 

It has been noticed (see, for example, |[T4l l20l ) that this sequence splits the 
natural numbers staying in the second line of Table |2]in two groups satisfying ([T|). 
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Table 2: The Morse-Hedlund sequence 



100101100 1 1 1 1 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 



Table 3: The half-completed magic square based on the Morse-Hedlund sequence 





2 


3 




5 






8 


9 






12 




14 


15 





The first group contains the numbers corresponding to ' 1 ' in the Morse-Hedlund 
sequence and the second group contains the numbers corresponding to '0'. Thus, 
for the first 16 numbers we obtain 

li+4i+6i+7i + 10i + ll' + 13i + 16i = 2I+3I+5I + 8I+9I + I2I + I4I + I5I, 
l2+42_,_62_,_72_^ ^02+112+132+ = 22 + 32 + 52 + 32 + 92+122+142+152, 
13+43 _^63_^ 73 _^ ^0^+113 + 133 + 163 = 2^V3^ + 5^ + 8^ + 9^'' + 123 + 143+15^ 

It can be seen immediately that these identities are a particular case of our results 
with 

a = l, b = 4, c = 2, d = 3, ki=0, = 2^ ytg = 2l (26) 

Naturally, by taking = 2" it is possible to generate Prouhet identities for any n. 

Let us now establish relations of the introduced method to the theory of magic 
squares (and cubes) developed in [T.Tl. AT xT matrix is called a magic square of 
order T if the sum of the entries in any row, column, or diagonal is the same num- 
ber. Usually, the entries of the magic square are taken from the numbers 1 , 2, . . . , 
and each of the numbers is used exactly once as an entry. 

There exist different ways to construct magic squares (see f3l). Let us consider 
the one proposed in li^J and based on the Morse-Hedlund sequence. In this paper 
we consider only the core case T = 4 being a basis for construction of the squares 
(and cubes) of the order T = 2" (see flllll for details). 

The algorithm to make the 4x4 magic square is the following. Count through 
the boxes in the square, and whenever the number you count corresponds to '0' in 
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Table 4: The completed magic square based on the Morse-Hedlund sequence 



16 


2 


3 


13 


5 


11 


10 


8 


9 


7 


6 


12 


4 


14 


15 


1 



Table 5: This matrix is a magic square for all real numbers a,b,c, and d satisfying 
dll) and any real numbers ki and k2 





a 


b 


c + ki+k2 


c + ki 


b + k2 


a + k2 


d + ki 


C + k2 


b + ki 


a + k\ 


d + k2 


d 


a + ki+k2 


b + ky+k2 


c 
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Table 6: The magic square with irrational entries obtained by using the data (|27] ) 



5 + V2 + V3 





6 


l + V2+\/3 


1+V2 


6 + n/3 


V3 




1+V3 


6 + V5 




5 + \/3 


5 


V2 + V3 


6 + V5+V3 


1 



the Morse-Hedlund sequence (see Table |2l), write the number down in the box. The 
rest of the boxes is filled in by counting numbers from 16 to 1 and gives the magic 
square presented in Table |4] 

Since the Prouhet identities were obtained as a particular case of our method 
with the data from (l26l ) used in construction of the magic square in Table |4l let us 
execute the inverse operation and substitute in Table |4] the numbers 1,2, . . . , 16 by 
the corresponding explicit formulae for the sets {jc, : 1 < / < 8} and {y, : 1 < / < 8}. 
The result is shown in Table |5] were, of course, the numbers a,b,c, and d satisfy 
©. 

It can be shown that, due to the sum of the entries in any row, column, or 
diagonal is equal to 2{a + b + ki + k2). Thus, the obtained matrix is a magic square 
for all real numbers a,b,c, and d satisfying ([5]) and any real numbers ki and ^2- Of 
course, some combination of these parameters can give repeating numbers. Note 
that the data from (l26l ) give us the the magic square from Table |4] as a particular 
case of Table [5] 

We conclude the paper by an example with irrational numbers. 

a = 0, b = 6, c=l, d = 5, ki = V2, k2 = V3. (27) 
For these data the corresponding systems ([T]i and (|4j become 

+ 5>+v^'+(6+v^)l + \/3'+(6+v^)'+(l + v'2+v'3)'+(5 + ^+^3)' = O' + 6' + (1 + ^2)' + (5 + v^)' + (I + v^)' + (5 + v^)' + (v^+ + (6+ v^) 
+ 5^ + v^- + (6 + ^/2f + (6 + V3)- + {l + V2+ s/3f + (5 + ^2 + v^f = 0^ + 6- + (I + s/lf + (5 + V2f + (1 + v^)^ + (5 + v^)^ + {^/2+ V?)^ + (6 + ^3) 

+ 5^ + (6 + \/2)3 + \/3V (6 + \/3)3+ (1 + ^2+^)3 + (5 + ^2 + 71)3 = 0^ + 6^ + (1 + v/2)3 + (5 + + (1 + + (5 + + (\/2+ ^3)3 + (6+ v^) 

and 

|-+52 + v/2- + (6 + v^)2 = 0- + 6- + (l + y2)2 + (5 + v^)2, 
+ 53 + V2V(6+v/2)3 + v/3V(6+v^)3+(l + V2+V3)3 + (5 + V2 + v^)3 = 0^ + 6^ + (1 + v/2)3 + (5 + ^2)^ + (1 + 73)^ + (5 + ^3)^ + (\/2+ ^3)3 + (6+ v5+ ^3) 

Finally, the magic square obtained from the data (|27] ) by using the general method 
from Table |5]is shown in Table (6] 
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4 Appendix 



This Appendix contains an alternative description of the sets {xf : I < i <k} and 
{yi : 1 < ' < ^} proposed by the unknown referee. 

Given n, we construct ^c, , 1 < / < 2"+^ , and yi,l <i < 2"+^ , so that for 1 < m < 
n + l, 

I = 1 (28) 

i=l i=l 

Choose a, ^, c, and J so that a + b = c + d. Choose any n values ki , . . . ,k„. The 
Xi are obtained by adding a to an even subset of the kj (including the empty set) or 
b to an even subset of the kj or c to an odd subset of the kj or d to an odd subset of 
the kj. The yi are obtained by adding either a or to an odd subset of the kj and c 
or d to an even subset of the kj. 



Equation (1281) is then equivalent to the statement that the function of a and b 
given by 



fia,b)= £ (-1)1^1 
5C{1,2,....«} 



iG5 



ieS 



(29) 



is a function of a + b when 1 < m < « + 1 . 
A simple expansion shows that 



fia,b) 



5 S^at ;=0 \ J J \ieS 



m-j 



m—\ 
m—\ 

l^{ai + bi) 

m— 1 

7=0 
m— 1 

l^{ai + bi) 



m 
j 

m 
j 

m 
j 

m 
j 



m-j 



5^0 



I (-1)1^1 I 
5={i|,...,i|5|}^0 !i+...+!|i-|=m-; 



m- I 



I 

!l+...+i„=m-i 



/I + ... + (,; = m-; 



m — ] 



\s\ . 

m 

7=1 



41 



m-j 



n^j(-ir 

;=i 



(-l)"(?i+l)! xr--xnia + b + i:}^iki), ifm = «+l, 
(— l)"(?i)! ^1 • • •/:„, ifm = w, 

0, if m < n. 
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discussions and the unknown referee for the alternative description of the sets {x,- : 
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